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ABSTRACT 


The  effect  of  an  anisotropic  reverberation  pattern  on 
the  performance  of  various  signal  processing  systems  for  linear 
sonar  arrays  was  studied.  It  was  determined  that  the  use  of 
clipping  at  each  array  element  could  seriously  affect  the  per- 
formance of  the  system  in  comparison  to  a similar  analog  system. 
A comparison  between  square  law  and  multiplier-correlator  de- 
tection schemes  indicated  no  significant  performance  difference. 


HISTORICAL  BACKGROUND 


In  order  to  utilize  the  planar  array  sonar  system  to 
its  fullest  capacity  it  is  necessary  to  determine  the  "best" 
choice  of  signal  processing  and  beam  forming  circuitry.  This 
decision  v/ill  be  ultimately  based  on  a combination  of  cost, 
feasibility,  and  performance  considerations.  To  this  end  we 
have  analyzed  the  performance  of  several  systems  which  can  meet 
the  engineering  feasibility  and  cost  requirements. 

The  studies  involved  in  this  problem  involve  two  areas 
of  research.  First,  what  sort  of  environment  should  we  consider, 
i.e.,  what  is  the  nature  of  the  signal  and  noise  that  will  be 
encountered  in  practice?  Second,  how  effective  is  the  circuitry 
in  extracting  the  signal  from  the  noise?  Before  proceeding,  let 
us  give  some  historical  background  in  these  areas. 

In  the  early  history  of  noise  analysis  (e.g. , Rice  (1)), 
random  noise  was  considered  to  be  a stationary  Gaussian  process. 

As  long  as  one  considered  only  single  receiver  systems,  this 
model  proved  to  be  adequate  for  practical  analysis.  Later  a theory 
of  multiple  receiver  systems  was  developed  (e.g.  Faran  & Hills  (2) 
and  (3)) , and  it  became  necessary  to  develop  a model  to  include 
noise  cross -correlation  effects,  which  required  some  picture  of 
the  distribution  in  space  of  noise  sources.  To  obtain  this  cross- 
correlation, Jacobsen  (4)  examined  homogeneously  distributed 
sources  which  led  to  an  isotropic  distribution  of  noise  at  the  re- 
ceiver system.  In  many  realistic  situations  the  isotropic  assump- 
tion breaks  down  and  non-isotropic  fields  must  be  considered 
(Bryn  (5)).  Furthermore,  it  has  been  found  in  practice  that  cer- 
tain non-Gaussian  effects  (sometimes  erroneously  called  non-stati^n- 
ary)  must  be  considered.  These  effects  include  random  amplitude 
modulation  (Thomas  and  Williams  (6))  and  random  phase  fluctuations 
(Berman  and  Berman  (7)). 


Paralleling  this  development  is  the  evolution  of  noise 
models  based  on  a more  physical  approach  to  the  problems  under 
consideration.  In  sonar  work  these  problems  appear  to  fall  into 
three  main  categories:  impulse  noise,  covering  a wide  range  of 
phenomena,  such  as  "static"  in  the  electronics,  wave  slaps, 
shrimp  in  Chesapeake  Bay,  etc.,  flow  noise,  and  reverberation. 
Impulse  noise  has  been  studied  by  Middleton  (8)  and  others,  prin- 
cipally in  relation  to  radar  and  communications,  rather  than  sonar. 
Flow  noise,  important  at  high  speeds,  is  related  to  pressure  fluc- 
tuations of  the  boundary  layer,  which  has  been  studied  by  Kraichnan 
(9),  Corcos  (10),  Gardner  (11),  and  others,  and  to  cavitation  at 
the  boundary.  Reverberation  has  been  examined  by  many  authors 
(e.g. , Stewart  (12),  Westerfield  (13),  Faure  (14)),  and  seems  to 
be  the  principal  problem  at  low  ship  speeds. 

When  we  examine  the  circuitry  we  can  for  the  purposes 
of  analysis  immediately  divide  the  elements  into  two  categories, 
linear  and  non-linear.  Linear  problems  involve  determining  the 
optimum  frequency  response  of  filters  (Childers  & Reed  (15)) , the 
number  of  elements,  geometry,  and  amplitude  and  phase  responses  of 
receivers  (Bryn  (5)),  methods  of  digital  sampling  (Rudnick  (16)), 
etc.  Non-linear  problems  concern  the  use  of  clipping  (Gore  (17) , 
Rudnick  (16)),  and  detection  instrumentation,  e.g.,  choice  between 
square  law  or  some  sort  of  mulitplier  correlation  scheme  (Faran  & 
Hills  (3)). 

In  this  paper  we  have  studied  the  non-linear  processing 
with  the  particular  emphasis  on  realistic  noise  models.  The  prob- 
lems of  particular  interest  are  the  effect  of  clipping  when  the 
noise  background  is  non-isotropic,  and  also  a comparison  under 
these  conditions  of  various  combinations  of  clipping  (non-clipping) 
and  detection  schemes. 
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SECTION  I 
INTRODUCTION 

In  this  note  we  are  concerned  with  the  problem  of 
determining  as  precisely  as  possible  the  angular  distribution 
of  noise  and  signal  with  the  objective  of  determining  if  sig- 
nals are  present  and  if  so  in  what  direction. 

Once  the  filter  of  the  detection  system  is  matched  to 
that  of  the  anticipated  signal,  the  only  characteristic  which 
distinguishes  signal  from  noise  is  the  fact  that  the  signal  is 
concentrated  in  one  direction,  so  that  it  would  appear  as  a local 
spike  in  the  angular  pattern  of  signal  plus  noise.  Furthermore, 
if  (as  in  the  case  of  an  active  pulse  system)  range  information 
is  also  available,  the  signal  can  be  seen  as  a spike  in  the  range 
distribution  in  a given  direction. 

For  this  problem,  a linear  array  of  hydrophones  is  used 
as  the  reception  system,  and  we  are  concerned  with  ascertaining 
the  efficiency  of  various  processing  systems.  In  particular  we 
are  concerned  with  (1)  the  effect  of  clipping  and  (2)  a comparison 
of  square  law  with  multiplier-correlator  detection,  for  an  arbi- 
trary (anisotropic)  noise  background. 

In  order  to  compare  the  various  systems,  we  have 
adopted  as  a figure  of  merit  the  signal/noise  ratio  (as  defined 
by  Faran  & Hills  (2))  i.e.,  the  ratio  of  the  difference  of  the 
on-target  and  off-target  mean  outputs  to  the  r.m. s.  value  of  the 
on-target  output.  Although,  in  general,  this  criterion  can  be 
very  misleading  in  many  situations,  e.g.,  if  the  smoothing  times 
were  different,  for  these  systems  (particularly  for  arrays  of 
many  elements)  we  can  use  the  results  for  qualitative  comparisons. 
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SECTION  II 


SYSTEMS  DESCRIPTIONS 


The  basic  hardware  in  the  system  under  study  is  a linear 

array  of  M (0 M-l)  hydrophone  receivers,  ^ach  followed  by  a 

band  pass  filter,  with  uniform  spacing  between  adjacent  elements. 

In  order  to  determine  the  direction  of  a signal,  each  array  element 
uses  an  associated  delay  device,  and  in  all  cases  analyzed  we  will 
assume  a uniform  spacing  of  the  delay  times,  i.e.,  the  n—  element 
(n  = 0,...,  M-l)  will  have  its  input  delayed  nA,  where  A is  some 
time  interval. 

Given  this  basic  setup,  we  are  concerned  with  comparing 
four  different  signal  processing  systems,  as  follows: 

1.  The  output  (i.e.,  the  delayed  signals)  are  added,  then 
the  resultant  is  squared  and  finally  smoothed  by  a 
simple  time  averaging. 

2.  The  array  is  divided  into  an  upper  and  a lower  half, 
(elements  (0,  . ..,  L-l) , (L,  ...  , M-l)).  The  out- 
puts of  each  half  are  summed,  then  the  two  sums  are 
multiplied  together  and  the  resultant  is  smoothed. 

3.  Same  as  1. , except  that  before  adding,  each  output  is 
infinitely  clipped,  i.e.,  only  the  polarity  informa- 
tion is  retained  at  each  hydrophone.  (Clipping  and 
delay  commute,  so  that  the  question  of  which  to  do 
first  is  determined  purely  by  hardware  considerations.) 
Furthermore,  since  clipping  generates  higher  harmonics, 
it  is  necessary  to  filter  the  clipped  signal  before 
squaring.  (See  Gore  (17)). 

4.  Same  as  3.,  except  that  after  the  second  filtering, 
the  array  is  divided  into  two  halves  and  the  final 
processing  is  the  same  as  in  2. 


SECTION  III 


NOISE  DESCRIPTION 

For  this  problem  we  assume  that  the  noise  sources  are 
at  distances  sufficiently  large  compared  to  the  dimensions  of  the 
array,  so  that  the  sound  waves  from  any  direction  are  parallel 
as  they  encounter  the  array.  However,  we  do  not  assume  that  the 
noise  is  isotropic  at  any  given  time.  Although  it  is  not  necessary, 
we  assume  the  array  is  backed  by  an  infinite  baffle  so  that  only 
sources  in  a hemisphere  are  received. 

Let  0 be  the  angle  between  the  direction  of  a noise 
source  and  the  array,  and  let  0 denote  the  angle  with  a vertical 
axis  (the  array  is  assumed  horizontal) , then  u(0,  0,  t)  will 
denote  the  pressure  at  time  t at  the  element  of  the  array 
arising  from  noise  in  the  (0,0)  direction.  Let  d denote  the  dis- 
tance between  adjacent  array  elements  and  c the  speed  of  sound, 
then  the  pressure  at  the  n—  element  is  then  u(0,  0,  t + — cos  0). 

Let  a (t)  be  the  total  sound  pressure  at  time  t at  the 
n—  element,  then 

an(c)  = 27  / / u(0,0,t  + ^ cos  0)  sin  0 d0  d0 


To  complete  the  model  we  will  assume 

E(u(0,0,t)u(0'  ,0'  ,t' ))  = 2irP(0 ,0) 6 (cos0-cos0' ) 6 (0-0' ) p (t-t ' ) 

where  P(0,0)  is  the  noise  power  in  the  (9,0)  direction  and  p(s) 
is  the  noise  autocorrelation  function  (which  we  are  assuming 
independent  of  direction  since  it  will  be  determined  by  receiver 
filter  characteristics) . The  anisotropy  in  the  noise  power  can 
be  considered  to  be  a product  of  anisotropy  in  source  and  in  trans- 
mission. 

When  a signal  is  present,  it  is  considered  as  noise, 
which  is  distinguishable  from  the  background  as  a delta  function 
in  some  direction. 


_ 
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SECTION  IV 


GENERAL  DERIVATION 


.th 


Let  an(t)  be  the  amplitude  of  the  pressure  at  the 
receiver  at  time  t,  let  be  the  output  of  the  k—  system,  then 
we  have  (where  A is  the  relative  delay  between  adjacent  receivers) 
k=l,  no  clipping -square  law 


G1  - T 


/ 


M-l 

2 

o 


k=2,  no  clipping  - multiplier 

T ;L-1 


1 1 /L_i 

g2  = t / S an(t“nA) 

o l o 


an^t_nA)l  dt 


M-l 


2 an(t-nA)  I dt 
L 


Let 


p(t,z(t))  = / f (t-x)  z (t)  d-r 
o 

be  the  output  of  a filter,  which  has  an  impulse  response  f(t) , with 
input  z(t).  Then  for 


k = 3,  clipping  - square  law 


= T f 


M-l 

p|t,  2 sgn 
o 


-12 

lan(t-nA)JJ  dt 


k=4,  clipping  - multiplier 


M-l 

an(t-nA)||  p(t,  2 sgn 


an(t-nA)||dt 


Let  us  assume  that  there  is  a signal  in  the  direction  0 . then  if 

d t h ^ 

AQ  = — cos  ©o,  the  mean  on- target  response  of  the  k—  system  is 

E(Gjc(Aq)).  To  obtain  the  off-target  response  there  are  two  alterna- 
tive approaches,  (1)  choose  a range  of  angles  near  0Q , sufficiently 
far  away  so  as  not  to  be  in  the  main  beam  steered  to  0Q,  and  take 
some  weighted  average  of  the  mean  responses  at  these  angles,  or  (2) 
assume  the  signal  is  absent  and  calculate  the  resultant  ECG^CA^). 


Although  the  former  is  a more  realistic  approach  for  most  practical 
situations,  the  latter  has  the  advantage  of  computational  ease,  and 
further  the  difference  between  the  results  of  the  two  procedures 
would  in  general  be  negligible.  Therefore  our  figure  of  merit 
is  given  by 


\ - 


E<Ck(V>  - E<Gk<V> 


where  indicates  the  output  when  the  signal  is  absent.  In  pass- 
ing, note  that  Faran  and  Hills  used  as  their  criterion  r£.  The 
immediate  problem  then  is  to  determine  ECG^)  and  aCG^)  as  functions 
of  A.  Let  us  introduce  some  additional  notations: 

Let  vn1,n2<tl't2>  “ E(ani(tl’an2(t2)> 

~n1,n2(tl-t2>  ‘ E|s«n[an1(tl)an2<t2>]l 


Since  we  are  dealing  with  Gaussian  variables,  we  have  (see  Appendix  B) 

^ 2-1  Tni ’n?(tl’t2)\ 

vV  (t,  ,t9)  = — sin  I — - — 1 

nl,n2  1 2'  tt  l Z I 

2 

where  Z = E(an(t)  ),  independent  of  n and  t,  is  the  mean  power  at 
any  element.  To  calculate  aCG^) , we  need  the  fourth  order  moments 


E an1<tl)an2<t2>an3(t3)an^V  = ^.n/h'V  V,  .n^h*  V 


(b-i  > t9)  v 


2 *4 


(t9«tA)  + v 


n1*n3'  1 * S'7  n9  ,nA  v 2 * 4'  T % ,nA  ^1  ’ vn9,n9^2’u3 


1 ’ 4 


(tl.tA)  VT 


2 >3 


(bo » to)  , 


(since  the  process  is  Gaussian).  For  the  clipped  case  in  general 

E jsgn  cannot  be  represented  exactly  (see  Placket  (18)). 

However'^in  ^he  large  array  case  since  we  are  concerned  only  with 
the  products  of  sums  of  the  form  Z a (t,  ) , where  the  number  of 

Hi  K 

nk  k 

terms  is  large,  we  may  assume  an  approximate  Gaussian  behavior  (see 
Tze-Chien  Sun  (19)) , at  least  to  the  extent  that  the  fourth  order 
moments  may  be  expressed  in  terms  of  the  second  moments,  i.e. , 


Let 


z*=  K a*n\\\\ 


Then 


E JT  zj**  E^Z^  E(Z3Z4)  + ECZ^g)  E(Z2Z4)  + ECZ^)  E(Z2Z3) 


o - 1 / > n 4 > n > 

E(zkZj)  - I 2 5 sin  H z ' J 

nk  nj  \ 


Then  for  k = 1,2 


T Ik“1  M-l 


E(Gk)  = Y / 2 Z vn>m(t-nA,  t-mA)dt 

? T T Ik'1  M-l  Ik"1  M-l 

0(0  - T"  / ; Z zz  z 

T o o o Jk  o Jk 

[vn1>n3<trnla-t2-n3a)  vn2,n4<trn2a't2-'V> 

* 

W*  vn2,n3<tl-n2a-t2-n3a)J  dtl  dt: 


where  I1  = M,  I9  = L 


Jk  " M-Xk 


For  k = 3,4 


T t t 

ECGfc)  - |r  / / / f(t-i2) 

O O O 


k M-l  l(v  (t,  -n/  ,T0-ia'  \ 

X Z sin-1  u'.  1 - ’ -2  ..L  dr.  dr,  dt 
o J.  z | 1 l 


»2  T T t2  C2  C1  fcl 


°2(Gk>  “ (It)  / / / / / / f(t1-T1)f(t1-T2)f(t2-T3)f(t2-Ti 

I o o o o o o 


dTldT2dT3dT4dtldt2 


where  I^I^,  l2=I^>  and  J^M-I^  (as  above).  Using  the  represents 
tion  given  in  Section  III 

vn,m(tl’t2)  = 7F  f f p(e>0)p  + cosej  sine  d0  d0 


Let  us  assume  that  the  received  spectrum  is  given  by 


p (t)  = / F(co)  cos  2tt  o-'t  dco 
o 


DI,J^Xo,t-)  = 7 L f F(a')  Hj  ^ j|co(x-xo)  ,a.'t|  dcodx 

o 00  )w  00  (X3  1 1 2j+l 

1,J  ° j=o  J 1 o o -1-1  J J l = 1 £ 


QOO 


rr~  d54  dcoi 


n=o  m=o 


cos  2t r v+  — (n-m+J)u 


where 


Q(x)  = f / P(e,0)d0 

o 


XJ 

j 2j+1  \ , 

J \i=i  i=j+i  I * * ° 

I 3 2j+l  \ 

fi.  — I 2 2 to. 

J lg=l  i=J+l  1 


Then  for  k = 1 or  2 


E(Gk>  ” DIk,\(xo'0) 


«V  - + “lk.Jk<*o.*>  DIk,Jk(x0,-t)] 


and  for  k = 3,4 


E<°k>  - EIk,Jk(Xo>0) 


o2<°k>  • 7 l <T-«  (E?k,0<xo-«  * EIk,Jk(xo-t>EIk,Jk<xo,-t>)  dt 

(after  expanding  sin  ^ in  a power  series  and  using  the  results 
in  Appendix  C) . 

We  can  simplify  the  expression  for  Ej,j(xo,t)  as  follows 


g (x  t)  = / / 

n °»  o -1 


1 e-27rico(n|(x-xo)+tj  p(flo)  9£gl  dxdo 


Then 


2 fl-1  1-1  “ l,2  , ,2j 

El,J<xo,t>  * > [„E0  £ JTT  I Wr<V'l 


ft  l~1  I_1 

7_  £ 

7 r n=0  m=0 


* B K-n^J^o,^  Re  Sn-m+J^o,^ 


SECTION  V 


ANALYSIS 


First  let  us  compare  the  two  cases  without  clipping. 
For  these  cases,  there  is  no  restriction  on  M.  We  observe  that 
E«V  is  a linear  functional  of  Q(x) , the  angular  distribution 
of  the  signal  plus  noise.  As  a result  we  see  that 

E(Gfc)  - E^*)  = EtGfc**)  » where 

JLA 

Gu  is  the  output  when  there  is  only  signal,  i.e., 


Q(x)  = 2S6(x-xq),  when  the  signal  is  at  xq 

Therefore  ECC^)  - ECG^*)  = l£  S , since  Hj  j (o,o)  = l£ 

k’  1c 

To  obtain  o^(Gx),  let  Q(x)  = 2S6(x-xq)  + QN(x)  , when  QN(x)  is 
the  angular  distribution  of  noise  alone,  then 

DI>J(xo,t)=sl2p(t>  + 7 / Qn(x)  / FC^Hj^  j^(x-XQ)  ,a)tj  daidx 


For  large  values  of  1,  the  noise  term  in  the  variance  is  propor- 
tional to  I and  the  noise  amplitude  in  the  direction  xQ,  so  that 
it  is  in  general  small  compared  to  the  signal  term,  therefore 


<Gk>^ 


“V 


(T-t)p*(t)dt 


so  that 


Rk* 


fy/T<T-t)p2(C)dt 


i.e.,  there  is  no  significant  difference  in  performance  between 
the  square  law  and  multiplier  correlator  schemes.  If  we  include 
the  noise  term,  then  is  increased  by  a small  term  proportional 
to  1/Ik»  which  implies  a slight  superiority  for  the  square  law 
system. 


For  small  M,  e.g. , M = 2,  the  noise  term  dominates  Dj  j 
(i.e. , for  k = 1,  1=2,  for  k = 2,  I = 1).  For  k = 1,  we  have 

1 °°  i d \ 

D2,o^xo,1:^  = Qn^x)  / F(“)(cos  2iru)t)  jl  + cos  2ww  - (x-xo)j  da)dx 


For  k = 2 we  have 


1 1 “ 

Dl,o(xo,t:)  = 7 f1  QN(x)  f F<“)  cos  27rQt  dudx  = NP(t) 

Dl  l(xo,t)  = j f Qn(x)  / F(w)  cos  2Tra)|t  + | (x-xQ)  dcudx 


We  can  simplify  by  observing  that  in  general 


i 1 fs  in]  , / , 1 \ 

j f Qn(x)  < or\  2tto£  (x-xQ)dx  « N = j f QN(x)dx 
-1  cos  I -1  J 


As  a result 


a2(G1)  = / (T-t)  p2  (t)  dt 

o2(C2)  = ^ / (T-t)p2(t)dt 
T o 


Therefore 


fl/l  (T-t)  (t)  dt 


R2  = 


Y~\/2  / (T-t)  p2  (t)  dt 
v o 


In  passing,  note  that  for  M = 2,  if  we  had  used  the  on  target 
versus  off  target  criterion,  would  assume  1/2  of  the  above 
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value.  Next  if  we  compare  k = 3 to  k = 4 for  large  M,  we  would 
get  the  same  result,  i.e.,  there  is  not  significant  difference 
between  square  law  and  multiplier-correlator. 


Finally  let  us  compare  the  unclipped  and  clipped  cases 
for  large  M,  i.e. , k = 1 and  3 compared  or  k = 2 and  4 compared. 
Essentially  the  difference  lies  in  the  difference  between  the  ex- 
pressions for  Dj.  j and  Ej  j.  The  key  difference  is  the  presence 
of  the  term  b||  8n_m+j(x0> t)||  *n  the  expression  for  Ej  j,  where  the 
corresponding 'term  in  j is  unity.  Therefore  let  us  summarize 
the  important  properties  of  B.  These  are  B(o)  = B(l)  = 1,  B(u) 
is  an  increasing  function  of  u for  0 £ u ^ 1 (B  is  defined  only  on 
the  unit  interval).  Next  let  us  examine  |sn(x0»t)|  . We  see  that 
this  is  the  magnitude  of  a weighted  average  of  complex  numbers  of 
unit  magnitude.  In  particular,  jgQ(x0»t)|  5 1.  If  Q(x)  is  con- 
centrated in  a given  direction  xq,  then  [sn(x0 !•  if  F(oo)  is 
very  narrow  band,  then  |gn(x0,t)|  is  independent  of  xq  and  t.  If 

Q(x)  is  constant,  then  g (x  . t)  = 0 for  n 4 0 (for  the  narrow 
cod  1 no 

band  case  when  — = •y) . To  see  what  this  means,  let  us  represent 
EI,J  h? 


EI,J**o,t)  ■4DI,J(xo't>  + -§T  /.  <3<x>  1 f(<i>)Gj  J(x0.t;x,o)da)dx 

ttZ-Io 


where 


I-l  I-l  / / \ I . , 

Glj(xo,t;x,a))=nZo  JE  B | gn_m+j(x0»0|J-  J cos  2ttu>  |t+(n-ra+J)-(x-xo)j 


First  let  us  examine  the  case  when  Q(x)  is  almost  uniform.  Then 

Gj  q«^I(1  - J)cos  27T<jot 

and 

EI,0<xo-t>*'  k + “ » - 

IT 

while 
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EI,I^xo»t^  T DI,l(xo,t:) 

V 

Therefore 

n “m.O^o’0^”  Hr  DM.0(xo’0^ 

R3^  rf  ■ ~ ~ - - : =2= 

T (T_t)(DM,0(xo’t)+M(S+N)4  " 1Mt>)  dt 

r r DL.L<xo>°>  " ¥ Dt.L<xo’°> 

tV2^T (T't)[K.0(V)+M<S+N)<^  ' X)  P (t))2-DL (L (x„ , t) DL _ ^x^-C)]  dt 

Using  the  approximations  to  j made  previously,  and 
assuming  S « N,  we  have 


1 + m <7  " !) 


1 + 1 + 


R4  R2  ‘ 


N 4 

We  see  then  that  if  = -j— g-(—  “ smflH  compared 

to  unity,  clipping  has  no  appreciable  effect.  However,  if  it  is 
large  then  there  would  be  a serious  degradation.  If  the  noise 
is  isotropic  then  Y^  would  be  small.  However,  in  cases  where  the 
medium  behaves  in  such  a way  that  the  signal  and  noise  powers  in 
the  given  direction  are  greatly  reduced  relative  to  other  directions 
(i.e.,  Q(xo)/lkS  is  still  small),  the  unclipped  systems  could  still 
detect  the  signals,  while  the  performance  of  the  clipped  system 
would  be  degraded. 
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To  illustrate  the  effect  of  clipping,  the  mean  outputs 
of  two  systems  (without  and  with  clipping)  using  square  law  detec- 
tion, were  calculated  for  a given  anisotropic  constant  frequency 
background  with  signals.  The  array  had  200  elements  spaced  .4 
wavelengths  apart.  As  seen  from  the  results  (Figure  1),  the 
apparent  noise  in  the  high  loss  direction  is  considerably  increased 
by  the  clipping. 


SECTION  VI 


CONCLUSION 


1)  Square  Law  Versus  Multiplier  Correlator 

The  analysis  indicates  that  from  a signal  noise  ratio 
point  of  view,  the  square  law  is  somewhat  better  than  the  multi- 
plier-correlator. However,  since  the  large  arrays  the  difference 
is  negligible,  it  would  be  necessary  to  make  an  analysis  based 
on  the  probability  distributions  of  the  output  to  determine  if 
either  system  has  any  significant  advantage. 

2)  Use  of  Clipping 

When  clipping  is  used  we  see  that  in  certain  circum- 
stances, e.g.,  when  there  is  a high  medium  loss  in  the  signal 
direction,  there  could  be  a serious  degradation  in  detectability. 
To  obtain  a more  quantitative  picture  of  the  situation,  the  out- 
put (square  law  detector)  power  distribution,  for  a particular 
input  power  distribution,  was  calculated  for  both  the  unclipped 
and  clipped  cases.  The  results  of  the  example  illustrate  the 
clipping  loss  quite  well.  (See  Figure  1.) 

Finally,  we  can  observe  that  the  non-linear  effect  of 
clipping  on  the  noise  is  similar  to  the  effect  described  by  Ogg 
(21)  on  the  interference  on  a weak  signal  by  a strong  signal  when 
clipping  is  used. 
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LIZED  (TOTAL  POWER  CONSTANT) 


COSINE 


APPENDIX  A 


FIGURE  OF  MERIT 


As  noted  in  the  introduction,  we  are  using  as  a figure  of 
merit  the  deflection  criterion  of  Lawson  and  Uhlenbeck  (22).  It 
must  be  observed,  as  Marcum  (23)  has  shown,  that  this  is  not  a 
good  quantitative  criterion  for  threshold  detection.  A more  ac- 
curate procedure  would  involve  first  selecting  a probability  of 
false  alarm,  then  calculate  a threshold  setting  from  the  system 
output  distribution  with  signal  absent,  and  finally  calculating 
the  probability  that  the  output  will  exceed  this  threshold  when 
the  signal  is  present.  In  practice,  this  procedure  is  extremely 
difficult  to  carry  out,  since  the  system  output  distributions  are 
difficult  to  determine  except  in  certain  special  cases. 

As  a result  we  are  forced  in  practice  to  use  the  deflec- 
tion criterion,  at  least  for  a qualitative  analysis,  i.e. , when 
the  signal  to  noise  ratio  is  small,  ^1  or  less,  then  systems  with 
significantly  different  ratios  will  perform  as  predicted  by  the 
criterion.  Also,  if  in  one  system  the  ratio  is  » 1,  and  in  an- 
other «1  or  less,  then  the  former  can  be  considered  superior. 
However,  if  two  different  systems  give  large  signal-to-noise 
ratios,  then  there  is  no  way  of  comparing  then  unless  something 
further  is  known  about  the  underlying  distributions. 


17 


APPENDIX  B 

CORRELATION  OF  CLIPPED  GAUSSIAN  VARIABLES 


Let  X^,  X2  be  correlated  Gaussian  variables  of  mean 
zero,  and  p = correlation.  We  wish  to  obtain  E(sgn(X-^ ^2) ) . 

Lemma  1:  E(sgn(U))  = 2 Prob(U  >0)  -1,  if  Prob  (U  = 0)  =0 

Proof:  Obvious 

Without  loss  of  generality,  we  may  assume  all  deviations  equal  one. 

2 - 1 

Theorem:  E(sgn(X^X2))  = — sin  (p) 

Prob(X1X2  2 0)  = 2 Prob  (Xx  0,  X2  2 °) 


- 00  00 

—F=?t ! e 

,Vi7 0 0 


2 2 

X1  + x2  “ 2pXj^ 
2(1-P2) 


dx^dx2 


it  _ r (l-2p  cos©  sinQ) 

1 00  "2"  TTt  27 

/ / e ( _P  ) r drd8 

-\/i  ? o o 

tt  v 1-p 


li£l  r7 


J l-2p  cos©  sin© 


= 7 cos-1(-p) 


1 , sin' 

7 + v 


Therefore,  by  Lemma  1,  the  theorem  is  proved. 
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APPENDIX  C 
FILTER  EFFECT 


Let  f(t)  = Coe~^a|t|cos  amn^t  be  the  impulse  of  the 


band  pass  filter  and  let 


2 j+1 

[j(t)  * E 


cos  2co  Trt  + (n,-n0)Pm 
m I L m 


We  wish  to  evaluate  Yj(t^,t2)  given  by 


C1  t2 


Yi(ti,t2)  = / / f(t1-T1)f(t2-T2)X.(T1-T2)dT1dT2 

-*  o o J 

By  a change  of  variables 

tl  t2 

Y.(ti,t2)  = / j f (s^) f (so)^^ (s2~sl  + t^-t2)ds^ds2 
j o o J 

Let  us  ignore  the  transient  effect  of  the  filter,  so  that 

00  00 

Y . ( 1 1 , t2)  - f f f (s^)  f (s2)^  ( s 2 — s i + t^-t2)ds^ds2 
-J  o o J 

00  00 

= / / f (s) f (s+u)X . (u+t^-t2) duds 

rt  _c  -J 


O OO  00  00 


= / / + / / f (s) f (s+u)X. (u+t^-t2)dsdu 

_ - 1 1 n n *■» 


-oo  -u  o o 


= / X.(u+t^-t2)  / f(s)f(s+  u|)dsdu 

— no  J n 


=>  c lf  X.  (u+t1-t2)  f(u)du 
— 00  J 


To  carry  out  the  integration  we  must  first  digress  slightly. 
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n 


: I 


\ 

s ; 


Let 


2n 

Z2n  * IT  cos  Ak 


i 2n 

( 


2n  / iAk  -iAk 
e + e 


02n  i Z e i A, 
i 2 k m,k  k 

r?n  mil  e 


* €m,2n^ 


where  e , = + 1,  such  that  the  vector  7 = (e 

m,k  — ’ 'm  m,l, 

takes  on,  as  a function  of  m,  all  possible  values.  Therefore 

02n-l 


J2n  ~ ;Tn-T  _Z,  COS (A1  + g €m»kAk 


m=l 


2n 


(For  simplicity  we  have  defined  7 so  that  € , = 1 for  m = 1, 

ry  1 in  m j 1 

2zn”  and  -1  for  the  remaining  m)  . 


Now  let  us  examine  a typical  term  in  Yj , using  the  Z2n 
representation;  such  a term  W is  of  the  form 


W = C / e~2na  U cos (Au+B) du 


= 2C  cosB  / c“2Trau  cos  Audu 


47raC  cos  B 

2 ;"T2 


(4ttcx)  + A 


Now 


2j  tl  / 

2ttgj  u j]"  cos  2a>ml  Tru+7r(t^-t2)  + (n^-^)^ 
m=l 


cos 


m 


,2j  + l 


Z cos 


n=l 


2ttu  (cjl> 


2j+l 


Z € _ CD 


o , n,m  m 

m=  i 


2 j+1 

+ 2 Z e cjj 

, n ,m  m 
m=l  ’ 


7r(trt2)  (nrn2)em) 


* 
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